CHAPITRE 6 



CARACTERISTIQUES STATIQUES DES PALIERS CIRCULAIRES 

EN REGIME LAMINAIRE 



1 - INTRODUCTION 

Les paliers lisses sont frêquemment utilisês ; les plus simples sont constituês d'un 
arbre cylindrique tournant a l’intêrieur un alêsage en prêsence de lubrifiant. Sur certains 
mêcanismes, ils correspondent a la meilleure solution technologique existant actuellement. On 
les rencontre dans les moteurs thermiques, les compresseurs, les turbomachines, les 
alternateurs, les rêducteurs, etc... La figure 1 prêsente diffêrents coussinets de paliers 
circulaires. 




Fig. 1 : Coussinets (Doc. LE COUSSINET MODERNE) 

Un palier lisse est constituê de deux êlêments, l'arbre gênêralement en acier de rayon 
R a et le coussinet en bronze ou en acier rêgulê de rayon intêrieur R c et de longueur L. Ainsi, 
un palier lisse peut être schêmatisê par deux cercles de rayons voisins et caractêrisê par trois 
grandeurs : le jeu radial C = Rc - R a , le rapport C/R (ou jeu relatif de l’ordre de 10' 3 ) et le 
rapport L/D de la longueur au diamètre du palier. La figure 2 prêsente les trois phases que l’on 
observe lors de la mise en route d’un palier sous charge. Les points 0 a et 0 C reprêsentent 
respectivement les centres de l’arbre et du coussinet et W dêsigne la charge extêrieure 
appliquêe a l’arbre. 

A l’arrêt, (fig. 2a), l'arbre et le coussinet sont en contact, la gênêratrice commune est 
opposêe a la charge et la distance 0 c 0 a est êgale au jeu radial C. 

Au dêmarrage, fig. 2.b, l’arbre roule en glissant a l’intêrieur du coussinet, le rêgime de 
lubrification est dit mixte, le lubrifiant est entrainê dans l'espace convergent formê par l'arbre 
et le coussinet. Dès que la vitesse de rotation devient suffisante, il y a crêation d'un champ de 
pression hydrodynamique qui s’oppose a la charge (fig. 2c). Pour une vitesse de rotation 



stable, et une charge W constante, le centre de l’arbre 0 a occupe une position fixe a l’intêrieur 
du coussinet. Si le torseur des forces extêrieures se rêduit a une force unique W qui agit dans 
le plan de la section mêdiane du palier, les axes de l'arbre et du coussinet sont parallèles. 
Deux paramètres sont alors nêcessaires pour positionner 0 a par rapport a 0 C . On utilise 
habituellement la distance 0 c 0 a appelêe excentricitê e et l'angle que fait la direction de la 

charge W avec la ligne des centres 0 c 0 a appelê angle de calage 4>. 






Fig. 2 : Schêmatisation d’un palier lisse 



Dans ce chapitre, nous nous limiterons au rêgime hydrodynamique êtabli et nous 
supposerons que l’êcoulement du fluide a travers le palier est isotherme, laminaire et 
pennanent. Le processus de calcul des performances d’un palier peut être divisê en 4 parties : 

- Ecriture de l’èquation donnant l’êpaisseur du filin lubrifiant, 

- Ecriture de l’èquation de Reynolds correspondant au palier êtudiê, 

- Dêfinition des conditions aux limites sur la pression, 

- Rêsolution des êquations et calcul des performances. 

2 - EOUATIONS DE BASE 

2.1- Epaisseur du film lubrifiant 



Considêrons une section droite du palier (fig 3a). Un point M appartenant a la surface 
du coussinet sera repêrê par la coordonnêe angulaire e-fcA , O c mJ ; l’êpaisseur du film 
d’huile est donnêe par : 

h = 0 C M - 0 C M’ = R c -O c M'=R a +C-O c M’ 



En appliquant la règle des sinus au triangle 0 a M’0 c il vient : 

0 a M’ e R a R a R a 

sin(0 a 0 c M') sin a sin(M'0 a 0 c ) sin(jt-0) sin 0 

d'ou 

sina = -^-sin0 et 0 C M' = -^— sin (o c 0 a M') 
R a sin 0 
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or 




Fig. 3a : Section droite du palier Fig. 3b : Palier dêveloppê 



en dêveloppant le sinus il vient : 



0 C M'= R... 1- sin 0 -ecosO 

c \ R. 



la quantitê e/R a < C/R est très petite devant l'unitê (de l’ordre de 10" 3 ) on peut donc nêgliger le 

/ e \ 2 

tenne sin 0 devant 1 . L'êpaisseur du li 1 m s'êcrit alors : 



h = C (l 



+ e cos i 



ou l’excentricitê relative : e = e/C varie de 0 a 1. II faut remarquer que cette relation n’est 
valable que si C/R « 1 . 

2.2 - Equations de Reynolds 



Les rayons R a et R c êtant très voisins, on peut nêgliger la courbure du film, dêvelopper le 
coussinet (ou l’arbre) et l'assimiler a un plan (fig. 3b). Les conditions cinêmatiques relatives 
aux surfaces s'êcrivent ; 

U, = R 0) ; V =W =0 

C C 11 

U =R oo cosa ; V =R oo sina ; W =0 

2 a a 2 a a 2 
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ou a, angle entre O c M’ et O a M’, est la pente de la courbe h(x) : 



ordre : 



dh dh e . 

tga = — = = sin 0 

dx Rd0 R 

la quantitê e/R a êtant très petite, on peut effectuer un dêveloppement et se limiter au T 

cosa «1 et sina ® a » dh/dx 



dans ces conditions les vitesses des surfaces s’êcrit : 



U, = R oo 

1 c c 

U = R co ; V = R oo — 

2 a a 2 a a dx 

Si on reporte les vitesses des surfaces du contact dans l’êquation de Reynolds êcrite en 
coordonnêe cartêsiennes et si on pose R = R a = R c , il vient : 




Remarques : 

a) On a dêveloppê le coussinet et choisi ce dernier comine plan de rêfêrence ; un rêsultat 
identique aurait êtê obtenu en dêveloppant l’arbre. 

b) Les vitesses angulaires oo a et oo c s'ajoutent. ainsi la portance est proportionnelle a la soinme 
des vitesses. Si celles-ci sont êgales et opposêes, la portance de l’êcoulement est nulle au 
second ordre près. 

2.3 - CONDITIONS AUX LIMITES 



Pour calculer le champ de pression, il faut rêsoudre l’êquation de Reynolds sur le domaine 
reprêsentê fig. 3b, compte tenu des conditions aux limites sur la pression : 



- selon z : 

- selon 0 : 



0, z = ^r) =P, ; 



L\ 



P| d ’ z = ~J 



=P. 



p (0 = 0, z^= p (0 = 2jt, z) Vz 



V0 



ouPi et P? reprêsentent les pressions extêrieures aux extrêmitês du palier. si le palier possède 
des rainures d’alimentation en lubrifiant, il faut alors êcrire les conditions aux limites 
correspondantes : 

pour une rainure circonfêrentielle en z = 0 : p (0, z = 0) = P a 
pour une rainure axiale situêe a l’abscisse 0 = 0 a : p(0 a , a) = P a 
pour une alvêole rectangulaire 0 (0 a , 0b) et z (z a , z c ) : p = P a 
ou P a est la pression dans la rainure ou dans l’alvêole. 

De plus, coinme la fonne du film est convergeante puis divergente en fonction de 0, il 
peut y avoir rupture du filrn dans le divergent lorsque la pression devient infêrieure a la 




pression de vapeur saturante du fluide. Les conditions aux limites relatives a cette situation 
ont êtê prêsentêes de facon gênêrale, dans le chapitre 3. 

Pour un palier hydrodynamique en fonctionnement stationnaire, trois types de 
conditions aux limites ont êtê successivement dclinies : 

2.3. a - Conditions de Sommerfeld 

Sommerfeld en 1904 [1] suppose que le film lubrifiant est continu et qu’il n’y a pas de 
rupture du filrn. Dans ces conditions la distribution de pression est antisymêtrique par rapport 
au point 0 = jt et le lieu du centre de l’arbre dans le coussinet est une droite normale a la 
direction de charge. 

Ces conditions aux limites ne sont valables que pour les paliers fonctionnant a de très 
faibles charges ou avec des pressions d’alimentation très êlevêes car sinon des pressions 
nêgatives apparaissent dans le filrn, ce qui est physiquement inacceptable. 

2.3. b - Conditions de Giimbel 

Gumbel en 1921 [2], a la suite des travaux de Sommerfeld, propose coinme conditions 
aux limites : 



p(0 = o, z) = 0 
p(0 = jt, z) = 0 
p(0, z) = 0 si jt<0<2jt 

Cela revient a nêgliger, dans la solution de Sommerfeld les pressions nêgatives. Ces 
conditions impliquent une discontinuitê dans l’êcoulement au point 0 = jt ; elles ne sont pas 
acceptables physiquement rnais sont utilisêes dans le cas de l’approximation du palier court. 

2.3.c - Conditions dites de Reynolds 

Ces conditions proposêes indêpendamment par Swift en 1932 [3] et par Steiber en 
1933 [4], respectent la continuitê du dêbit et supposent que pour une abscisse 0 S inconnue, la 
pression et le gradient de pression s’annulent : 



p(0 = 0, z) = 0 

P(0 = 0 S , z) = 0 

^(0-e„z)-iP<e-0,.z)-o 

00 d z 

p(0, z) = 0 si 0 s < 0 < 2jt 

On distingue ainsi deux rêgions dans le palier : une rêgion active ou la pression est 
positive et le filin complet, pour des valeurs de 0 comprises entre 0 et 0 S et une rêgion 
inactive a pression nulle et dans laquelle en rêgirne dynamique, peut apparaitre de la 
cavitation, pour 0 s < 0 < 2jt. Ces conditions sont celles qui donnent les solutions les plus 
exactes ; elles ont êtê vêrifiêes dans de nombreuses êtudes expêrimentales. 




2.4 - Rêsolution de l'êquation de Reynolds 



L'êquation de Reynolds est une êquation aux dêrivêes partielles du second ordre de type 
elliptique qui n’admet pas, sauf dans quelques cas simples, de solutions analytiques. 

Sommerfeld [1] a rêsolu cette êquation dans le cadre d’un problème plan (5/5 z = 0) en 
nêgligeant l’êcoulement axial, (hypothèse du "palier long"). Une autre solution 
unidimensionnelle, moins satisfaisante du point de vue analytique mais plus reprêsentative du 
point de vue physique a êtê proposêe par Ocvirk [4] dans le cas ou la longueur du palier est 
très faible devant son diamètre ("palier court"). 

Pour les paliers de longueur finie, diffêrentes mêthodes de calcul ont êtê proposêes : analogie 
êlectrique, mêthodes variationnelles, dêveloppements en sêrie et mêthodes numêriques. Les 
mêthodes actuellement utilisêes sont pour la plupart des mêthodes numêriques aux diffêrences 
(inies et aux êlêments (inis. 



3 - LE "PALIER INFINIMENT LONG" 



Lorsque le rapport L/D de la longueur au diamètre du palier est grand on peut nêgliger 
le dêbit axial. L'êcoulement est donc purement circonfêrentiel. En pratique, ce calcul 
s'applique aux paliers tels que L/D > 4. Si seul l’arbre est animê d’un mouvement de rotation 
uniforme oo a = w et oo c = 0, l’êquation de Reynolds (5) se rêduit a : 



_d_ 

dx 



h 3 



dp N 

dx y 



= 6 q Ro) 



dh 

dx 



ou en utilisant la coordonnêe angulaire 0 



x 

R 



de [ de ) 



6 p R 2 (u 



dh 

dŌ 



3.1 - SOLUTION DE SOMMERFELD 



Supposons le palier alimentê a la pression p = P a par l’intermêdiaire d’une rainure 
axiale de largeur nulle (infiniment mince) situêe a l’abscisse 0 = 0. Les conditions aux limites, 
connues sous le nom de conditions de Sommerfeld s'êcrivent : 



p (0 = 0) = p (0 = 2jt) = P a 



Une première intêgration de l’êquation de Reynolds donne : 



dp 

d0 



6 p R 2 oo 



h-h : 

h 3 
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h* êtant l’êpaisseur du film lubrifiant au point ou le gradient de pression est nul. En intêgrant 
et en rcmplacant h par sa valeur, on obtient : 



R 



p = 6 [x tn — JJ 



d0 



h ! 



C) J (l + scos0) C J (l + ecos0) 






de 



. + K 



ou K est une constante d’intêgration. 



Pour calculer ces deux intêgrales, on peut soit utiliser les tables prêsentêes par Booker 
[5], et donnêes en annexe I, soit utiliser le changement de variable de Sommerfeld dêfini par : 

l-£ 2 



1 + SCOS0 = 



1 - e cosrjt 



soit 



cos 0 = 



COSlf - £ 
1 - £ COS 



sin 0 = 



(l-r) 



/2 



sinif 



1 - e cos 



d0 = 



(l — £ 2 ) drjt 

1 - s cos rjt 



Dans ce cas, la pression s'êcrit : 



ôlicoIR/C) 2 \ . h* 2xb -4 £ sinrlt + £ 2 xlt + £ 2 sinrlt cosilt ] 

P= 7 V v../ J x|t - £ sin x|t f \\ -l + K 

(l-e 2 ) 1 2 Ct-e 2 ) j 



Les deux constantes h et K sont dêterminêes a partir des conditions aux limites sur la 
pression. II vient : 

h* 2 (l - e 2 ) 



K = P a ; — = 



C 2 + C 



soit : 



P = 



6 p, (o (R / c) 2 
(l-e 2 )' 1 



hf - e sin - 



2xf - 4 e sinif + e 2 + e 2 sinif cosip 
2 + e 2 



+ p. 



La figure 4a reprêsente schêmatiquement l’allure de la courbe de pression. Celle-ci est 
antisymêtrique par rapport au point ip = n (ou 0 = Jt). On remarque que, selon la valeur de la 
pression d’alimentation P a , il peut y avoir ou non une rêgion ou la pression est nêgative. 

Pour calculer la charge portante W et l’angle de calage (p, il faut intêgrer le 
champ de pression sur la surface de l’arbre. Pour cela, on projette les forces de pression selon 
deux directions orthogonales k et t dont l’une est colinêaire a la ligne des centres (fig. 4b). 

Si on nêglige les contributions des contraintes tangentielles, les êquations d’êquilibre 
s’êcrivent : 



?jr 

W cos (f> + Lj o p cos 0 Rd0 = 0 
- W sin (p + L f p sin 0 Rd0 = 0 
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Fig. 4a : Allure de courbe de pression 
En intêgrant par parties : 




Wcosd) = LRf — sin0d0 
d0 

W sin 4» = LR f — cos 0 d0 
d0 



pour calculer ces intêgrales, on effectue le changement de variable donnê prêcêdemment. II 
vient : 



d'ou : 



W cos cj) = 0 



W sin (|) 



W = 12 p. 



12 p, R 3 o) L it £ 

C 2 (2 + e 2 )Vl-s 2 

R 3 L 71 £ 




II est d'usage en lubrilication d’introduire le nombre sans dimension S appelê nombre 
de Sommerfeld : 



q LRoo / R 
K W (c 

La charge W, est par consêquent, le nombre de Sommerfeld S, sont indêpendants de la 
pression d’alimentation P a car celle-ci intervenant connne une constante dans la valeur de la 
pression admet une rêsultante nulle après intêgration autour de l'arbre. La ligure 5 prêsente les 
variations de la charge sans dimension W en fonction de l'excentricitê relative £. L'angle de 



\ 2 _ Vl-e 2 (2 + b 2 ) 
, 12 jr 2 £ 



1 19 



calage 4> est toujours êgal a ji/2. Cela signifie que le centre de l’arbre se dêplace selon une 
direction normale a la ligne de charge. On peut montrer qu'un tel palier est toujours 
dynamiquement instable. 




Fig.5a : Variation de la charge sans 
dimension W en fonction de l’excentricitê 
relative £ 




Fig.5b : Variation de l’angle de calage 4» 
en fonction de l'excentricitê relative £ 



Fig. 5 : Influence des conditions aux limites 



Le couple de frottement est obtenu par intêgration des contraintes de cisaillement a la 
surface de l’arbre. 



c "- r!l /o" T »f- h ) de 



avec 



xy 




et 



1 dp 
2q dx 



y(y-h) + ^ 

h 



d'ou 



T 



xy 



Roo 1 dp 
h 2 dx 



(2y-h) 



Le couple s'êcrit : 



C a = R 2 L 




Rw h dp 1 
h 2 dx J 



d0 
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soit encore : 



_2n 



C = R 2 L J’jT ■! M. — - + — — — 



R(o C dp ecos0 dpl 

u + — + L d0 

h 2R d0 2R d0j 



La valeur de la première intêgrale est immêdiate a partir des tables donnêes par 
BOOKER [5]. La deuxième intêgrale est nulle, tandis qu'en intêgrant par parties on reconnait 

Q 

dans la troisième intêgrale, l’expression — W sin c() . Le couple de frottement sur l’arbre s'êcrit 
donc : 



c a = 



R 3 q oo L 2 jt 
C 7TT7 



+ — W sin ò 
2 



soit compte tenu de l’expression de la charge : 



c = 



4 q (D R 3 L Jt 
C 



1 + 2 e' 



(2 + 8 2 ) 



Les variations du couple sans dimension C a =C a /(_LDNR(R/C))en fonction de 
rexcentricitê relative 8 sont prêsentêes fig. 6c. Un calcul semblable au prêcêdent permet de 
dêtenniner le couple de frottement sur le coussinet. 




Fig. 6 : variations du couple sans 
dimension C a en fonction de l’excentricitê 
relative e : Influence des conditions aux 
limites. Cas du palier infiniment long 



c = 

soit 



R 3 [X o) L 2 
”” C 



— W sin d) 
2 



c = 



4jtq o) R 3 L 
C 




On remarque que : 

C - C = e W sin <\> 

Cette relation traduit l’êquilibre du 

palier. 

Le coefficient de frottement relatif 
a l’arbre est : 

f a = C /(RW) = — J— sin cj) + ■ 

R 2 V7I7 



Afin d’êliminer le rapport C/R, on utilise le nombre de frottement : 



f=C/(CW) 



R 

C^ 



, . A 2jt 2 S 1 + 2e 2 
= 2 sin 9 + . + 

VI - e 2 3e 
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3.2 - SOLUTION DE GUMBEL 



La solution de Sommerfeld n’est valable que pour des paliers très faiblement chargês 
ou alimentês sous forte pression car sinon on obtient des pressions nêgatives dans le film 
lubrifiant. A la suite des travaux de Sommerfeld, Giimbel propose de nêgliger les pressions 
nêgatives. Ceci introduit une discontinuitê du dêbit dans le film (chapitre 4). Bien que 
physiquement et mathêmatiquement inacceptable, cette solution est frêquemment employêe 
car elle permet d’obtenir dans le cas du "palier court" par exemple, une solution proche des 
rêsultats expêrimentaux et elle conserve a l’êquation de Reynolds son caractère linêaire. 

Dans le cas d’un "palier long" alimentê a l’abscisse 0 = 0 avec une pression 
d'alimentation êgale a la pression atmosphêrique, c’est-a-dire P a = 0 , il vient : 

W cos 4> = -12 q (o RL (R/ C 
Wsinc)) = 6pcoRL(R/C) 2 



1 (-«’)&♦«’) 

e jt 

V 1 - ( 2+e2 ) 



Soit encore : 



W = 



6 qooRL (R/C) 2 e 
(2 + e 2 ) (l-s 2 ) 



^4 8 2 + JT 




jrVl-s 2 

<t> = arctg — 

2 8 



Les variations de la charge W et de l’angle de calage 4> en fonction de l’excentricitê 
relative sont prêsentêes sur les figures 5a et 5b. Le calcul du couple rêsistant est effectuê en 
intêgrant les contraintes de cisaillement sur l'arbre. Cependant, l’existence d'une zone active 
remplie d’huile (0 < 0 <jt) et d’une zone inactive (jt < 0 < 2jt) partiellement vide oblige a 
dêcomposer le calcul du couple en deux parties. Pour cela, on suppose que dans la zone 
inactive, le dêbit circonfêrentiel reste constant et êgal a celui calculê a la sortie de l’arc actif, 
c'est-a-dire dans ce cas a l’abscisse 0 = n. Cette hypothèse permet de dêfinir une largeur 
èquivalente L' dclinie par la relation : 



Q c = 



R o) h (0 = jt)l 
2 



Ro)h(0)L' 

2 



ou 




si h s reprêsente l’êpaisseur du film a l’abscisse 0 = jt. 



II vient ainsi : 



C, = R 2 L 



p, R o) 
h 



h dp' 

+ — 

2 dx 



d0 + 



f 



2 ^ q R co h s 



h 2 
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soit après intêgration : 



et 



C a = 



C c = 



jt [i 00 R 1 L 


(2 + 8 3e 2 

_j_ 


cVl-e 2 1 


(1 + 8 2 + 8 


jt p 00 R 3 L j 


(2 + 8 


3e 2 


C Vl-e 2 


1 + 8 


2 + 8 



Les variations du couple de frottement sur l'arbre sont prêsentêes sur la figure 6. 
3.3 - SOLUTION DITE DE REYNOLDS 



Pour des conditions d’alimentation identique a celles utilisêes prêcêdemment, c'est-a- 
dire une pression d’alimentation nulle a l'abscisse 0 = 0 (ou r|i = 0), nous savons que, s'il y a 
rupture du film lubrifiant, la continuitê du dêbit impose la continuitê du gradient de pression 
le long de la frontière qui sêpare la zone active (p >0) de la zone inactive (p = 0) 

(chapitre 4). Dans le cas unidimensionnel, cette condition s'êcrit : 




(e=e s ) 



= 0 



ou 0 S est l’abscisse du dêbut de la rêgion inactive. Sa valeur est une inconnue du problème. 

Dans le cas envisagê, la frontière ^ ou le film se reconstitue est connue. Celle-ci 
est situêe sur la rainure d’alimentation. 



Les deux constantes d’intêgration qui apparaissent dans l’êquation de pression sont 
donc dêtenninêes a partir de deux conditions aux limites : 



P (0 = O) = 0 ou p (rp = 0) = 0 



p (e = ° s )=f%) = 0 ou p (ip = h» s ) = ( -7^ 



d0 fi= 



0=0, 



dxf 



= 0 






On obtient alors le champ de pression : 



pM = 



6[iR (o 
C 2 (l-e 2 ) 2 



- 8 sin - 



2xf - 4e sin + s 2 + s 1 sin cos 



2 (l - 8 cos s ) 



ou rp s , qui correspond a l’abscisse de rupture du film d’huile, est dêfini par la relation 
8 (sin ip s cos s -'ip s )+2(sinxf s -rp s cosrf s )=0 



par intêgration du champ de pression la charge s’êcrit : 




W = 3q RooL 



R\ 

C 



V 1 - e 2 (l - 8 cos xp s ) « l -8 



e 2 (l-cosxb Y . / . v 

2 + 4 (sin xp s - xp s cosipj 



l'angle de calage est donnê par : 



tg <t> = 



2 Vl — s 2 (sinrlJ s -ip s cosx[) s ) 
e (l - cos s y 



le nombre de Sommerfeld S s’êcrit alors : 



S = 



VI - e 2 (l-ecosif s ) 



„ e 2 (l - cos . / . v 

3jt a 2 ~ + 4 ( sm if s - if s cos rp s / 

V 1-e 



Le couple de frottement est obtenu par l’intêgration des contraintes de cisaillement sur la 
surface de l'arbre, ou du coussinet. Coinme prêcêdemment, le domaine est sêparê en une zone 
active (0 < ip < if s ) et une zone inactive (if s < if < 2n). 

Sur l’arbre le couple s’êcrit : 



c a = 



p, R 3 oo L 
C 



2n + e(sinif s -if s cosif s ) 



V 1 - £ 2 (l — E COS if s ) 



+ — W sin 4> 



sur le coussinet : 



c = 



p, R 2 oo L 
C 



2n + e(sinif s -if s cosif s ) 



V 1 - £ 2 (l-SCOSlf s ) 



— W sin ò 
2 



Le nombre de frottement relatif a l’arbre est donnê par la relation : 



f=C a /(CW) 



= n S 



2jt + s(sinif s - if s cosif s ) 

V 1 - S 2 (l-SCOSlf s ) 



sin4) 



3.4 - COMPARAISON 



Les rêsultats obtenus avec les conditions aux limites de Sommerfeld, de Giimbel et de 
Reynolds sont prêsentês sur la figure 6. On remarque que les solutions de Giimbel et de 
Reynolds sont voisines, l’êcart le plus faible est obtenu pour les grandes excentricitês, ce qui 
est tout a fait logique car dans la solution de Reynolds, l'abscisse de rupture du fdm lubrifiant 
tend vers n lorsque l’excentricitê relative 8 tend vers 1. La figure 6 montre que l’effet du 
gradient de pression sur le couple ne devient sensible que pour des excentricitês relatives 
supêrieures a 0,5. 
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3,5 - Cas d’une position quelconoue de la rainure d’alimentation 



Dans la situation envisagêe prêcêdemment, nous avons supposê que la rainure 
d'alimentation êtait situêe au point d’êpaisseur maximale du film lubrifiant. Cette hypothèse 
qui simplifie les calculs, n’est pas rêaliste car dans la pratique la rainure d’alimentation est fixe 
vis-a-vis de la direction de la charge (on ne s'intêresse qu'au cas stationnaire). 



Considêrons le cas d’une 
alimentation sous pression nulle par une 
rainure d’alimentation situêe a cp degrês en 
aval de la direction de la charge (fig. 7). 

Les solutions obtenues a partir des 
conditions aux limites de Sommerfeld ou 
de Giimbel ne sont physiquement pas 
acceptables car elles impliquent soit de 
l'existence de pressions nêgatives, soit une 
discontinuitê dans le dêbit. La solution 
correcte est celle qui fait intervenir les 



frontières de 



reconstitution 



(S 



rupture 
1 



fc' 



) 



et la 



J 



du filrn lubrifiant. 




Ce problème a êtê traitê de manière rigoureuse [6]. Nous ne rappellerons ici que les 
principaux rêsultats. On choisit l'origine de la coordonnêe circonfêrentielle a la rainure 
d’alimentation ; l’êpaisseur h du filin lubrifiant est alors donnêe par : 

h (0) = C [1+8 cos (0 + a)] 
ou 

a = cp- 4> ; 0 G (0, 2 jt) ; a G (0, 2jt) et r G (0, 1) 

A la rainure d’alimentation et pour une pression d’alimentation êgale a la pression 
atmosphêrique nous avons : 

p(0) = p(2jt) = 0 



Dans la rêgion active Q+ la pression est positive et est obtenue a partir de l’èquation de 
Reynolds : 



d 

d0 



h 3 



dp N 

dê y 



= 6 pi R 2 co 



dh 

d0 



Dans la rêgion inactive £2o la pression est nulle (p = 0), la conservation du dêbit 
massique impose : 

d 



d0 



[j3(0) h(0)]= 0 soit (3 h = K = Cste 



les conditions aux limites pour la frontière de rupture ^ + s'êcrivent : 



P = 



dp 

d0 



= 0 ; (3 = 1 
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Les conditions aux limites pour la frontière de rêformation ^ sont : 

P - 0; hRm(l-|3)---L— ^ 

6 [x R d0 

Un certain nombre de rêsultats ressortent de ces expressions [4, 5] : 

a) la rêgion active fU est toujours connexe Q+ = (0i, 0 S ) ; 

b) pour rêsoudre le problème (3 (0) doit être donnê en un point de Qo. La connaissance de 
(3 entraine celle du dêbit massique par unitê de longueur dans Qo : 
point 0 = 0; 

c) (3(0) ne peut être quelconque. Si (3(0) est trop important, 
respectêe ; si (3(0) est trop petit, la pression ne peut s'êtablir. On a : 



n k. o) , • • , 

p(3 . On cnoisit le 



la condition a) n’est pas 



~ ~ < (3(0) < 1 

1 + 8 cos a 



Compte tenu de ces remarques, le champ de pression est alors facilement obtenu. En 
effet, si (3(0) h(0) est donnê on en dêduit immêdiatement l’abscisse 0 S de rupture du film : 

h(0 s ) = (3(0) h(0) 

puis le champ de pression par une intêgration numêrique de : 



dp 2 h(0)-h(0 s ) 

— = 6 [i R o) 



d0 



h J (0) 



pour 0e(-oo,0 g ) 



avec p(0 s ) = 0 

L’abscisse 0 = 0) de rêformation du film est le premier point ou p(0) = 0. 



Les rêsultats numêriques montrent que si l’abscisse de rupture 0 S dêpend peu de la 
valeur de (3(0), il en est tout autrement de l’abscisse de rêfonnation 0; [6] ; cela signifie que, 
dans le cas d’une pression d’alimentation êgale a la pression ambiante, le dêbut de l'arc actif 
est conditionnê par la valeur du dêbit d’alimentation. L'emploi des conditions aux limites dites 
de Reynolds, associê a l'algorithme de Christopherson (chapitre 3) n’impose pas de se fixer le 
dêbit, celui-ci est un rêsultat de calcul. Physiquement, cela revient a supposer que l’on est 
capable de fournir au palier toute la quantitê de fluide nêcessaire. Une comparaison [7] des 
rêsultats obtenus avec ces deux types de conditions aux limites montre que la charge portante 
est peu modifiêe tant que la rainure d'alimentation est situêe dans la zone non chargêe du 
palier (grand jeu) a condition de prendre pour (3(0) la valeur maximale possible. 

Si la formulation du problème est simple pour le cas unidimensionnel, il en est tout 
autrement dans le cas bidimensionnel pour lequel, le problème n’est encore que partiellement 
traitê [8]. 
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4 - LE "PALIER COURT" 



Lorsque le rapport L/D de la longueur au diamètre du palier est faible, on peut nêgliger 
le gradient de pression circonfêrentiel devant le gradient de pression axial. Cette hypothèse 
êmise pour la première fois de Michell [9] a êtê reprise et dêveloppêe par Ocvirk et Dubois 
[4]. Elle est entièrement justifiêe pour les paliers dont le rapport L/D est infêrieur ou êgal a 
1/8. En pratique cette hypothèse est encore utilisêe pour des rapports L/D de 0,5 car l’erreur 
commise reste faible sur le couple et le dêbit, elle est d’un ordre de grandeur acceptable sur la 
charge mais très importante sur la valeur du maximum de pression dans le film. Ces erreurs 
diminuent lorsque l’excentricitê dêcroit. 

L'êquation de Reynolds s'êcrit : 



^/ h3 dp\ 

dz \ dz 



= 6 q Rcn — 
dx 



-lh>A 



soit encore : 

dz \ dz 

Avec les conditions aux limites 



\ , dh 

= 6iift) — 

j de 



p| 0 ’ z = -y] = o et p ( e ’ z = + ^ 



= o 



Le champ de pression s'êcrit : 



p(e,z)-L^ z>-- 



2 \ 



8 sin 0 



(i + £ cos e ) 3 



II est êvident que la pression varie comme sin0 ; elle est positive pour 0 G (0, jt). Pour 

obtenir la charge, Ocvirk et Dubois utilisent les conditions aux limites de Gumbel ; il vient : 

0 2 

'R\ £ 



W = u LRoo — 
D 



/L\ 



C 



(l-£ 2 ) 



Vl6 



2 2 

8 +Jt 



d-e") 



Le nombre de Sommerfeld est tel que : 



et l’angle de calage : 



s-|A 

L 



(l ~ ) 

Tt EyJ 16 £ 2 + Jt 2 (l - £ 2 ) 



tg <t> = 



Jt V l-£ 2 

4 £ 



Le couple de frottement calculê en tenant compte des zones active et inactive du 
palier, est donnê par la relation : 
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Le dêbit axial est : 



Q z =LCVs 

ce qui justifie l’emploi d'un dêbit sans dimension 

Q = Q/LCV = e 



On remarque que le nombre de Sommerfeld et le nombre de frottement sont fonction 
du rapport (L/D) 2 ; ainsi, si l’on veut reprêsenter sur un seul tableau les caractêristiques du 
palier, il est nêcessaire de considêrer les nombres modilics de Sommerfeld et de frottement : 



S* = S| — 

D 



f * - 




f 



Les courbes de la figure 8 mettent en êvidence les êcarts comrnis sur le nombre de 
Sommerfeld S c’est a dire sur la charge adimensionnêe, lorsque cette hypothèse est retenue. 

5 - LE PALIER DE LONGUEUR FINIE 



Dans le cas gênêral, la difficultê d’obtenir une solution satisfaisante pour un palier n’est 
pas seulement dus a la forme de l’êquation de Reynolds mais êgalement aux conditions aux 
limites car le lubrifiant n’est pas toujours admis par l’intermêdiaire d’une rainure axiale a 
pression ambiante. Considêrons un palier lisse alimentê par une rainure axiale pratiquêe sur la 
ligne des centres dans la section d’êpaisseur maximum du film. 

L’êquation de Reynolds s’êcrit : 



d | 


fh> «E] 


i d i 

+ — 


fh’^1 


dx 1 


9xJ 


' dz 


9zJ 



= 6 p Rw 



dh 

dx 



Cette êquation n’admet pas de solution analytique ; diverses mêthodes approchêes ont 
êtê proposêes. La plupart sont dêrivêes des approches infiniment court ou infiniment long [10 
a 13] et ne permettent pas d’appliquer les conditions aux limites de Reynolds. Leur utilisation 
est souvent limitêe aux arcs partiels. Rohde et Li [14] ont proposê une approche qui pennet de 
rêduire l’êquation de Reynolds a une êquation diffêrentielle du second ordre en x. Utilisant le 
fait que la rêpartition de pression dans la direction axiale est parabolique dans le cas du palier 
court et reste encore quasi-parabolique dans le cas du palier de longueur finie, il posent : 

r, /_,2- 

p(0, z) = f(e) 



L 



et reportent cette relation dans la fonctionnelle associêe a l’êquation de Reynolds. 
L’intêgration par rapport a z, puis l’application du calcul des variations, conduit a l’êquation 
diffêrentielle : 
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_d_ 

d0 



If df' 
6fx d0 



10 

~6~ 



R 2 h 3 f(0) = 5 co 



dh 

d0 



Cette êquation doit être intêgrêe numêriquement. II est possible d'utiliser les conditions 
aux limites de Reynolds. Cette approche donne d'excellents rêsultats dans le cas des paliers 
circulaires, ou des arcs partiels, en rêgimes permanent et même en rêgime transitoire. Les 
rêsultats prêsentês sur le tableau 1 concement le cas du palier circulaire en rêgime pennanent ; 
l’emeur effectuêe sur la portance est donnêe en %. 

Tableau 1 : Approche de Rhode et Li [12]. Cas du palier circulaire : erreur sur la portance (%) 



8 

L/D^^ 


0,5 


0,8 


0,95 


0,5 


~0 


~0 


~0 


1 


~0 


1 


5 


2 


1 


4 


6 



Si l'on souhaite une plus grande prêcision, ou si la forme du palier ne permet pas de 
faire l'hypothèse d’une pression parabolique. II est nêcessaire de faire appel aux mêthodes 
numêriques. La plus simple est la mêthode des diffêrences finies, la mêthode des êlêments 
linis et êgalement employêe. 



Si on reporte dans l’êquation de Reynolds les variables sans dimension suivantes : 
0 = x/R, Z = z/L, h = h/C et P = p/6 p, oa (r/C) 2 

II vient : 



d 

dQ 



d0 



/ 



f— V— 

[l) ôz 



h j ^i 



dZ) 



dh 

d0 



Cette êquation ne fait intervenir que deux paramètres : rexcentricitê relative 8 et le 
rapport R/L c'est a dire L/D. Dans le tableau 2 et sur les courbes des figures 9 a 12, nous 
avons prêsentê les variations du nombre de Sommerfeld, de l’angle de calage, du nombre de 
frottement, du dêbit axial, de la charge et du couple relatif a l’arbre en fonction de 
rexcentricitê relative 8 pour diffêrents rapports L/D. Ces rêsultats sont insuffisants pour 
dessiner un palier car il faut aussi faire intervenir les donnêes de la rêsistance des matêriaux, 
êtudier la stabilitê, prendre en considêration les effets thermiques etc.... Cependant, ils 
permettent de faire un avant projet et d’obtenir un ordre de grandeur pour le dêbit, les pertes 
par cisaillement et la tempêrature moyenne du lubrifiant. Ces rêsultats dêcoulent d'une 
rêsolution numêrique de l’êquation de Reynolds par la mêthode des diffêrences finies. Le 
système linêaire obtenu est inversê de manière itêrative en utilisant l’algorithme de Gauss- 
Seidel avec surrelaxation ; les conditions aux limites sont celle dites de Reynolds. 



no 





Tableau 2 : Caractêristiques des paliers lisses circulaires alimentês a pression atmosphêrique 
et dont la rainure axiale est situêe sur la ligne des centres 

Notations (rappel) 



c 


Jeu radial C = R, - Rj (m) 


W 


Charge appliquêe N 


D 


Diamètre du palier (m) 


£ 


Excentricitê relative £ = e/C 


L 


Longueur du palier (m) 


<l> 


Angle de calage 


N 


Vitesse de rotation de l’arbre (t/s) 


q 


Viscositê dynamique (Pa s) 


Q 


Dêbit axial (m / s) 


00 


Vitesse angulaire de l’arbre (rd/s) 


R 

V = 


Rayon du palier (m) 
ooR Vitesse linêaire de l’arbre (m/s) 


c a 


Couple de frottement sur l’arbre (mN) 




Couple de frottement adimensionnê sur l’arbre & d 


= C a C/(qLDNR 2 ) = i|f a = | 



f Nombre de frottement sur Tarbre f = — f a = C d /(CW) 



f a Coefficient de frottement sur Tarbre f a = C d /(RW) 

Q Dêbit axial sans dimension Q = Q/(LCV) 

S Nombre de Sommerfeld S = (—) = ^ LDN ( — 

Jtw [ c } w [ c J 



"Palier court' ' : L/D < 1/8 



8 


0,1 


0,2 


0,3 


0,4 


0,5 


0,6 


0,7 


0,8 


0,9 


0,95 




0,99 


0,461 


0,272 


0,17 


0,106 


0,0625 


0,033 


0,0139 


0,00331 


0,000812 




83 


75 


68 


61 


54 


47 


39 


31 


21 


15 


-d-] 

c (dJ 


18,75 


8,514 


4,98 


3,14 


2,016 


1,25 


0,722 


0,355 


0,114 


0,0388 


e d 


18,94 


18,47 


18,31 


18,5 


19,02 


20,02 


21,89 


25,55 


34,58 


47,79 



1 


L/D = 1/4 


1 


8 


0,1 


0,2 


0,3 


0,4 


0,5 


0,6 


0,7 


0,8 


0,9 


0,95 


s 


16,2 


7,57 


4,49 


2,83 


1,78 


1,07 


0,58 


0,263 


0,0728 


0,0221 




82,5 


75,5 


68,5 


61,5 


54 


47 


39,5 


31,5 


21,5 


15,5 


R 

V* 


307 


140 


82,5 


52,67 


34,26 


21,85 


13,18 


6,97 


2,70 


1,20 


Q /( LCV ) 


0,0983 


0,196 


0,295 


0,393 


0,491 


0,590 


0,688 


0,787 


0,885 


0,933 


ê d 


18,95 


18,49 


18,37 


18,61 


19,24 


20,42 


22,74 


26,5 


37,09 


54,3 



ni 





L/D = 1/2 



8 


0,1 


0,2 


0,3 


0,4 


0,5 


0,6 


s 


4,32 


2,03 


1,21 


0,784 


0,508 


0,3180 




82 


75 


68,5 


61,53 


55 


48 


R 

L' 


82,10 


37,71 


22,55 


14,75 


9,94 


6,67 


Q /( LCV ) 


0,0938 


0,187 


0,281 


0,374 


0,468 


0,562 




19 


18,57 


18,64 


18,81 


19,57 


20,97 



L/D = 1 



8 


0,1 


0,2 


0,3 


0,4 


0,5 


0,6 


s 


1,33 


0,631 


0,388 


0,260 


0,178 


0,120 




79,5 


74 


68 


62,5 


56,5 


50,5 


R 

L' 


25,36 


11,87 


7,35 


5,07 


3,67 


2,70 


Q /( LCV ) 


0,0801 


0,159 


0,237 


0,314 


0,390 


0,466 




19,06 


18,81 


18,94 


19,5 


20,62 


22,5 



L/D = 2 



8 


0,1 


0,2 


0,3 


0,4 


0,5 


0,6 


s 


0,559 


0,271 


0,173 


0,122 


0,0893 


0,0654 




75 


71 


67 


62,5 


58 


52,5 


R 

L' 


10,76 


5,21 


3,40 


2,50 


1,96 


1,60 


Q /( LCV ) 


0,0538 


0,104 


0,153 


0,199 


0,243 


0,285 




19,25 


19,22 


19,65 


20,49 


21,95 


24,46 
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"Palier long" : L/D >4 
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Fig. 1 1 : Dêbit axial en fonction de l'excentricitê 

excentricite relative c 




L/D= 1 

L/D = o o 



Fig.12 : Angle de calage 
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6 - EFFET THERMIQUE GLOBAL 



L'hypothèse d’une viscositê constante en tout point de l'êcoulement est erronêe car on sait 
que la tempêrature est variable. La rêpartition de tempêrature a l'intêrieur du filin peut être 
obtenue en rêsolvant simultanêment l’êquation de l’ênergie, l’êquation de Reynolds gênêralisêe 
et les êquations de conduction de la chaleur dans l’arbre et le coussinet. La rêsolution de ce 
problème tridimensionnel n’a êtê envisagêe qu'a partir des annêes 1960 et une solution 
satisfaisante a ce problème n'a êtê donnêe qu'a partir des annêes 1980 [15]. Historiquement, et 
pratiquement, on a introduit le concept de "viscositê effective" qui consiste a supposer le 
rêgime isotherme et la tempêrature êgale a Tf = T e + AT, ou T e est la tempêrature d’entrêe du 
lubriliant et AT l’accroissement de tempêrature obtenu en êgalant un pourcentage de l’ênergie 
gênêrêe par cisaillement a la quantitê de chaleur êvacuêe par le lubriliant : 

P = L C a co = p C P Q AT 

ou A est un paramètre fixê : 0,5 < /. < 1 qui rend compte de l’environnement du palier. Pour un 
êcoulement adiabatique >, = 1. On prend souvent A= 0,85. Cela signifie que 85% de la chaleur 
est êvacuêe par le fluide et 15% par l'arbre et le coussinet. Coinme, pour une charge donnêe, 
le couple et le dêbit sont fonction de la viscositê du lubrifiant, l’êquation ci dessus doit être 
rêsolue de manière itêrative. Le processus est prêcisê fig. 13. 




Fig. 13 : Calcul de la viscositê effective 



La convergence de ce processus itêratif dêpend de la valeur de T e , de la forme de la 
relation p(t), du couple et du dêbit. Lorsque le processus diverge, on peut introduire un 
coefficient de pondêration o (sous relaxation 0 < a < 1) et êcrire [16] : 

AT (n+1) = (1 - J AT (n) + _ (_ C a )/( C p Q) 
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En pratique, une valeur de a êgale a 0,5 et une valeur initiale nulle pour AT assurent une 
convergence rapide. Cette mêthode a êtê introduite dans un grand nombre de processus 
automatiques de calcul sur ordinateur. 

Cette approche ne permet pas de dêterminer la tempêrature maximale de surface du 
palier. Certains auteurs [17] proposent d'utiliser pour cette valeur la relation empirique : 

Tmax = T e ° + 2 AT 

Si la tempêrature maximale ainsi calculêe est proche de la tempêrature limite 
admissible, une analyse thennique plus prêcise du palier doit être effectuêe. 

7 - APPLICATION AU CALCUL D'UN PALIER LISSE 

Soit un palier lisse lubrifiê avec une huile dont la viscositê, la masse volumique et la 
capacitê caloriliquc a la tempêrature de fonctionnement, sont respectivement : 

viscositê q = 0,015 Pa.s. 

masse volumique p = 860 kg/m 3 

chaleur massique Cp = 2000 J/kg C° 

Les dimensions du palier sont : 

longueur L = 100 nnn 

rayon R = 50 min 

jeu radial C = 0,075 mm 

la charge appliquêe est W = 50 000 N. 
et la vitesse de rotation N = 50 tr/s. 



Pour ces conditions, le rapport L/D = 1 et le nombre de Sommerfeld est : 



pLDN / R 



W C 



— =0,0666 



la vitesse linêaire d’un point appartenant a la surface de l'arbre est donnêe par: 

V = R o) = 15,71 m/s 
Après interpolation dans le tableau 2, on obtient : 

- excentricitê relative : e = 0,733 

- angle de calage : 4> = 4 1 ,3° 

- nombre de frottement : f — = 1,794 

C 

- couple sans dimension : = 27,63 

- dêbit sans dimension : Q = 0,569 
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On en dêduit : 



- L'êpaisseur minimale du li 1 m : h m j n = C (1 - s) = 20 iim 

- le couple de frottement : C x = 6,9 m N 

- le dêbit axial : Q = 67 10' 6 m 3 /s. 

- la puissance dissipêe : P = C a w = 2140 W 

Si l’on suppose que 85% de cette puissance a êtê êvacuêe par le fluide (k = 0,85) 
l'êchauffement moyen de l'huile est de : 

AT = _P /(_C P Q) = 15,8 °C 

Cela signifie que la tempêrature d’entrêe du lubrifiant est d’environ 16° infêrieure a 
celle correspondant a la viscositê utilisêe pour les calculs. 

Les solutions analytiques donnêes par les approches du "palier court" et du "palier 
long", conduisent aux rêsultats dimensionnês suivants (tableau 3): 

Tableau 3 : Comparaison entre "palier court", "palier long" et palier de longueur finie 





L/D = 1 


"Palier court" 


"Palier long" 


excentricitê relative 


0,733 


0,59 


0,38 


angle de calage 


41,3° 


46,2 


63,6 


nombre de frottement 


1,794 


1,34 


1,54 


couple sans dimension 


27,63 


19,9 


21,8 


dêbit sans dimension 


0,569 


0,59 


- 



On constate que les erreurs peuvent être très importantes, particulièrement sur 
l’excentricitê relative, c'est a dire sur l’êpaisseur minimale du filrn et sur l’angle de calage qui 
sont fortement surestimês. Elles restent modêrêes pour le couple de frottement, et 
pratiquement nêgligeables sur le dêbit dans le cas du palier court. 

8 - C OMP ARAISON THEORIE - EXPERIMENTATION 

Les mêthodes de calcul isothermes ne peuvent être vêrifiêes qu'en utilisant de l’eau 
comine lubrifiant ou pour des conditions de fonctionnement pour lesquelles l’ênergie dissipêe 
par cisaillement est faible et conduit a un êchauffement et une diminution de viscositê peu 
importantes. 

Des essais ont êtê rêalisês a l'INSA de Lyon [18] avec des paliers circulaires de 
diamètre 50 mm, lubrifiês a l’eau. On obtient une bonne corrêlation entre les valeurs mesurêes 
et calculêes, fig. 14 a 15 

Des êcarts significatifs sont constatês même pour des paliers de faible diamètre 
lorsque les hypothèses de calcul ne sont plus vêrifiêes. Les mêthodes de calcul isothermes ne 




conduisent pas a une bonne corrêlation avec les rêsultats expêrimentaux. C'est le cas, par 
exemple, pour une lubrification a l'huile, des rêsultats obtenus par DUBOIS et OCVIRK [16] 
pour un palier circulaire de diamètre 35 nnn (fig. 16 et 17). 
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Fig. 14a : Palier long Fig. 14b : Palier court 

Fig. 14 : Variations du nombre de Sommerfeld avec l'excentricitê d’après [15] 
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Fig. 16 : Angle de calage d’après [16] 
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ANNEXE 1 



Calcul des intêgrales du type [5] : 



T<m _ 

n = J(1 



sin' 0 cos m 0 



+ e cos 



avec 0 < e < 1 



Ces intêgrales s’êcrivent : 



(n - 1 ) (1 - e 



— £ sin 0 



+ £ cos 



- + (2n - 3)1" - (n - 2)1“ n * 1 



T oo 2 1 - £ 0 

•i = , v Arctg L~ tg- 

J1-p 2 Vl + e 2 



T oo 1 [ -esin0 _ t00 

2 = l-£ 2 (l + £COS0) 1 



2(l- 



- £ sin 0 (4 + 3 £ cos 0 - £ 2 ) 



+ ECOSt 



(4£ 2 -l)l[ ( 



t 01 _ j_ r T oo joo 

1 n — L A n + A n-1 
£ 



T 02 _ r T 00 _ ry T 00 T 00 

1 n — 2 cn Z1 n-1 + A n- 



t °3 _ _j_ f_ T 00 t T°0 _ T 00 T 00 

1 n - 3 L A n + ^n-l J X n-2 + X n- 



(n - i)e (l + £ cos 0) 1 



Y ; pour n * 1 



l[° =- — tn (l + £ cos 0); pour n = 1 
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